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ABSTRACT 
> 

^^ \ Sky masking is unavoidable in wide- field weak lensing observations. We study how 

pg \ masks affect the measurement of statistics of matter distribution probed by weak grav- 

.^ I itational lensing. We first use 1000 Gaussian simulations in order to examine in detail 

^^ ■ the impact of mask regions on the weak lensing Minkowski Functionals (MFs). We 

consider actual sky masks used for a Subaru Suprime-Cam imaging survey. The masks 

increase the variance of the convergence field and thus the expected values of the MFs 

Lj . are biased even for a Gaussian random field. The bias is caused by two effects. One is 

;^ \ owing to the reduced number of sampling Fourier modes, which can be accounted for 

- - -' analytically by considering the survey geometry appropriately. The other is owing to 

variation of the variance of the convergence field for each field of view. Lensing MFs 

are biased systematically when the reconstructed convergence field is normalized by its 

variance. We then use a large number of cosmological ray-tracing simulations in order 

to address if the lensing MFs measured for a Subaru 2-deg^ survey are consistent with 

those of the standard cosmology. The resulting x^ /^doi = 29.6/30 for combined three 

MFs, obtained by taking the mask effects into account, suggests that the observational 

data are indeed consistent with the standard ACDM model. Finally, we explore how 

masked covariance matrices affect cosmological parameter estimation. We conclude 

that the lensing MFs are powerful probe of cosmology only if the effect of masking is 

correctly taken into account. 



INTRODUCTION 



An array of recent obs e rvations such as the cosmic microw ave background (CMB) anisotropics 



[e^E:. 



Komatsu et al. 



Tegmark et al. 



2006 



2011 



Reid et al. 



Planck Collaboration et alj 120131 ) and the large-scale structure (e.g. 



2010l ) established the standard ACDM model. The energy content 



of the present-day universe is dominated by dark energy and dark matter, and the primordial den- 
sity fluctuations, that seeded all rich structure we observe today, were generated through inflation 
in the very early universe. There still remain a few important questions such as the nature of 
dark energy, the physical properties of dark matter, and the exact mechanism that generates the 
primordial density fluctuations. 



Gravitational lensing is a powerful method to study matter distribution (e.g. lOguri et al.ll2012l ). 
Future weak lensing surveys are aimed at measuring cosmic shear over a wide area of more than 
a thousand square degrees. Such observational programmes include Subaru Hyper Suprime-Cam 
(HSC) ^, the Dark Energy Survey (DES) ^, and the Large Synoptic Survey Telescope (LSST) ^. 
Space missions such as Euclid and WFIRST are also promising. The large set of cosmic shear data 
will enable us to improve the constraints on cosmological parameters which provides important 
clues to the mysterious dark components. 

A variety of statistics are proposed to characterize the large-scale matter distribution. Minkowski 
Functionals (MFs) are among the most useful statistics to extract non-Gaussian informations from 
a two-dimensional or three-dimensional field. For example, MFs have been applied to the observed 
CMB maps and pr ov ided comparable c o nstrai nts t o those obtained using the CMB bispectrum 
(JHikage et al.l l2008l ) . iMatsubara fc JainI (l200in and ISato et al.l (|200ll ) studied Om-dependence of 
weak lensing MFs. More recently, iKratochvil et al.l (|2012l ) showed that th e lensing MFs contain 
significant cosmological information, beyond the power-spectrum, whereas IShirasaki et al.l (|2012l ) 
showed that weak lensing MFs can be used to constrain the statistical properties of the primordial 
density fluctuations. 

These previous studies on weak lensing MFs consider often idealized cases. There exist, 
however, many observational effects in real weak lensing measurement. For example, imperfect 
shape measurement due to seeing and optical distortion, selection effect s of galaxies, uiicertai n 
redshift distribution of galax ies due to photometric redshif t error fe.g. iBolzonella et al.M200d ). 
noise-rectification biases (e.g. iKaiseiJ l2000l : lErben et al.ll200ll : iHirata &: Seliakll2003l ). and compli- 
cated survey geometry due to mask regi ons. Some of these effects on cosinic she ar power spectrum 
analysis have been already studied (e.g. iHuterer et al.l 120061 : iHikage et al.ll201lh . A comprehensive 
study of observational effects on lensing MFs are also needed in order to fully exploit the data from 



^http://www.naoj.org/Projects/HSC/j _index.html 
^ ht t p : / /www . darkenergy survey, org/ 
^http://www. lsst.org/lsst/ 



upcoming wide cosmology surveys. 

In the present paper, we study the impact of mask regions on measurement of weak lensing MFs. 
Masking effect can be one of the major systematics because MFs are intrinsicahy morphological 
quantities. We use a large set of numerical simulations to critically examine the effect of masking. 
We then directly measure the lensing MFs from real observational data obtained from a Subaru 
survey. We compare the observed MFs with the results of our ray-tracing simulations that include 
the effect of mask regions explicitly. 

The rest of the present paper is organized as follows. In Section [21 we summarize the basics 
of MFs and how to estimate MFs from observed shear field. In Section [3l we describe the data 
used in this paper and the details of numerical simulations which includes ray-tracing simulations 
of gravitational lensing. In Section [U we show the results of the impact of mask regions on lensing 
MFs. We also represent the comparison with observed MFs and ray-tracing simulation results. 
We then perform the simple analysis to characterize the impact of mask regions on cosmological 
constraints from lensing MFs. Concluding remarks and discussions are given in Section [5l 



2. MINKOWSKI FUNCTIONALS 

2.1. Basics 

Minkowski Functionals (MFs) are morphological statistics for some smoothed random field 
above a certain threshold. In general, for a given D-dimensional smoothed field S , one can 
calculate D + 1 MFs Vi. On §^, one can thus define 2-1-1 MFs Vq, Vi and ¥2- For a given threshold, 
Vo, Vi and V2 describe the fraction of area, the total boundary length of contours, and the integral 
of the geodesic curvature K along the contours, respectively. MFs are defined, for threshold v, as 

Vo{u) ^ 1. f dS, (1) 

where Qy and dQy represent the excursion set and the boundary of the excursion set for a smoothed 
field u{0). They are given by 

Q, = {e\u{e)>v}, (4) 

dQy = {e\u{e) = v]. (5) 

For a two-di mensional Ga ussian random field, one can calculate the expectation values for MFs 



analytically (JTomitalll986l ): 
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2.2. Estimation of Lensing MFs from Cosmic Shear Data 

We summarize how to estimate lensing MFs from observed shear data. Let us first define the 
weak lensing mass maps which corresponds to the smoothed lensing convergence field k: 



/C(6>)= I d^cp K{e-cj))U{ct)), 



(9) 



where U is the filter function to be specified below. We can calculate the same quantity by 
smoothing the shear field 7 as 



]C{e) = ld^jt{ct>:e)Qt{cj>), 



(10) 



where 7^ is the tangential component of the shear at position cf) relative to the point 0. The filter 
function for the shear field Qt relates to U by 
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We consider that Qt is defined with a finite extent. In this case, one finds 



u{e) 



dO 



,Qt 



Qt{0), 



(11) 



(12) 



where Oq is the outer boundary of the filter function. 

In the following, we consider the truncated Gaussian filter (for U) as 
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for 9 < 60 and U = Qt = elsewhere. Throughout the present paper, we adopt 6g = 1' and 
60 = 15'. Note that this choice of 9g corresponds to an optii nal smoothing scale for detection of 
massive galaxy clusters using weak lensing with ^source = 1-0 (JHamana et al.ll2004l ). 



We follow iLim &: SimonI (|2012l ) to calculate the MFs from pixelated JC maps. We convert a 
weak lensing field /C to x = (/C — {IC))/ao where ctq is the standard deviation of a weak lensing field 



on a map. In binning the thresholds, we set Ax = 0.2 from x = —5 to x = 5. W e have checked 



that the binning is sufficient to reproduce the analytic MFs formula (JTomitalll986l ) for 1000 maps 
of Gaussian random fields on a 5° x 5° map without mask. When we estimate lensing MFs on a 
/C map with mask, we discard the pixels within 29g from the mask boundaries, because tC data on 
such regions are affected by the lack of shear data. 



DATA 



3.1. Suprime-Cam 

We use the i'-band data from the Subaru/Suprime-Cam data archive SMOKA^. The obser- 
vation is characterized as follows. The area is contiguous with at least four pointings, the exposure 
time for each pointing is longer than 1800 sec, and the seeing full width at h alf-maximum (FWHM ) 
is better than 0.65 arcsec. The data is dubbed "COSMOS" in Table Al in iHamana et al.l (|2012l l. 



We conservatively use the data only within 15 arcmin radius from the field center of Suprime- 
Cam, because the point spread function (PSF) becomes elongated significantly outside of the central 
area, which m ay make PS F correction inaccu rate. Then mosaic stac king is performed with S CAMP 
(JBertinl [2009 ) and SWarp (JBertin et al.ll2002l) We use SExtractor (JBertin fc ArnoutsI [1999 1 and 
hfindpeaks of IMCAT software (JKaiser et al.lll995l ). and then the two catalogues are merged by 
matching positions of the detected objects with a tolerance of 1 arcsec. 



For weak lensing anal ysis, we follow the KSB method (JKaiser et al.l Il995l : iLuppino &: Kaiser 



1997: Hoekstra et al, 



19981 ). Stars are selected in the standard way by identifying the appropriate 



branch in the magnitude half-light radius (r/i) plane, along with the detection significance cut 
S/N > 10. We found that the number density of stars is ~ 1 arcmin"^. We use the galaxy images 
that satisfy the following three conditions; (i) the detection significance of S/N > 3 and z^ > 10 
where i' is an estimate of the peak significance given by hfindpeaks, (ii) rh is larger than the stellar 
branch, and (iii) the AB magnitude is in the range of 22 < i' < 25 ( where MAG_AUTO g iven by 
the SExtractor is used for the magnitude and slightly different from lHamana et al.l (|2012l )). The 
resulting number density of galaxies rigai is then 15.8 arcmin"^. We measure the shapes of the 
objects by getshapes of IMCAT, and correct for the PSF by the KSB method. The rms of the 
galaxy ellipticities after the PSF correction is 0.314. 

Next, we define data and mask regions by using the observed positions of the source galaxies 
as follows. We map the observation area onto rectangular pixels of width 0.15 arcmin. For each 
pixel, we check if there is a galaxy within 9d = 0.4 arcmin from the pixel center. If there are 
no galaxies, then the pixel is marked as a mask pixel. After performing the procedure for all the 
pixels, the marked pixels are masked regions, whereas the other pixels are data regions. However, 
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Fig. 1. — The reconstructed lensing field /C from the Subaru Suprime-Cam data. The convergence 
/C is computed from the ehipticity of 102342 source galaxies by Eq. (jlOp . The masked survey area 
(black portion) covers 0.34 deg^. The grey-scale bar shows the value of (/C — (/C))/o"o. 



we unmask "isolated" masked pixels whose surrounding pixels are all data pixels. 

Weak lensing convergence field fC is computed from the galaxy ellipiticity data as in Eq. (jlOp 
on regular grids with a grid spacing of 0.15 arcmin. The resulting mass map includes mask regions 
as shown in Figured! The masked regions cover 0.34 deg in total. The unmasked regions are found 
to be 1.79 deg . Note that we use only 0.575 deg in unmasked regions for lensing MFs analysis 
because we remove the ill-defined pixels within 29q = 2 arcmin from the boundary of mask. 



3.2. Gaussian and Ray-tracing Simulation 

We first use 1000 Gaussian convergence maps to study the impact of mask regions on the lens- 



i ng M Fs. We generate the maps for our fiducial cosmology using the fitting formula of ISmith et al 



(J2003l ) to calculate the matter power spectrum P{k; z). We then integrate the matter power spec- 
trum over redshift z, convolved with a weighting function for the source redshift -Zsource = 1- Each 
convergence map is defined on 2048^ grid points with an angular grid size of 0.15 arcmin. 



In addition to the Gauss i an sim ulations, we use 1000 ray-tracing simulations of weak gravi- 
tational lensing of ISato et al.l ( 20091 ) Q|. The ray-tracing simulations are performed on light-cone 
outputs that are generated by arranging multiple simulation boxes. Briefly, small- and large- volume 
A^-body simulations are placed to cover a past light-cone of a hypothetical observer with an aii gular 
extent 5° x 5°, from redsh ift z = to z = 3.5, similarly to the methods in lWhite &: Hul (J2000l ) and 
Hamana fc Mellierl (l200ll ). We set the source redshift ^^source = 1 for the ray-tracing simulations. 
The grid size and hence the angular resolution are kept the sam e as our Gaussian simulations 
described above. Details of the ray-tracing simulations are found in lSato et al.l (|2009l ). 



It is well-known that the intrinsic ellipiticities of source galaxies induce noises to lensing shear 
maps. We model the noise by adding random ellipiticities drawn from a two-dimensional Gaussian 
to the simulated shear data. We set the root-mean-square of intrinsic ellipiticities to be 0.314 and 
the number of source galaxies is set 15.8 arcmin" . The values are obtained from the actual weak 
lensing observations described in Section 3.1. 



4. RESULT 

4.1. Masking Effect on MFs: Gaussian Fields 

We first discuss the effect of masking on the lensing MFs. We calculate the MFs for 1000 
Gaussian convergence maps with and without mask. Figure [2] compares the results for Vq, Vi, and 
V2- In each upper panel, the black points are the results with mask whereas the gray points are for 
those without mask. The error bars show the standard deviation of the MFs measured from our 
1000 maps. We also plot the difference between the masked maps and the unmasked maps in the 
lower panel. Clearly the expected values of the MFs are biased with the mask. The amplitude of the 
bias corresponds to approximately 5 percent differences of the variance amplitude (see Appendix 

ED. 

It is important to study the effects of mask regions in a more quantitative manner. The variance 
of the MFs in each bin is larger for the masked maps, because the number of pixels used to evaluate 
the MFs is effectively reduced by the mask. We can approximately evaluate the increase of variance 
by scaling the effective survey area. The bias of weak lensing MFs is actually caused by two distinct 
effects. One is owing to the reduced number of sampling Fourier modes in a /C field because of 
the mask. This changes the amplitude of MFs because both Uq = ((/C — (/C))^) and af = (|V/Cp) 
get smaller. We find that the differences of the MFs between the masked Gaussian maps and the 
unmasked ones are indeed explained by systematic decrease in (Tq and ai (see Appendix |A] for 
details). Another effect is caused by variation of the variance uo among different fields of view. 
In order to calculate the lensing MFs in real observations, we need to normalize the reconstructed 



^For the simulations, the adopted cosmology is consistent with WMAP3 year results i Spergel et al.ll2007l ) 
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Fig. 2. — Masking effect on the lensing MFs. In the upper panel, the black points with error bars 
show the MFs with mask whereas the grey points are the results without mask. The lower panels 
show the differences of MFs between the masked and unmasked maps. We use 1000 Gaussian 
convergence maps with Gaussian smoothing with a length scale of 1 arcmin. The error bars in the 
lower panel indicate the standard deviation estimated from the 1000 maps. 



lensing convergence field by its variance. The normalization then causes fluctuations of the MFs 
from expected values because the MFs are defined as a function of x = {IC — (/C))/cro where (Tq is 
evaluated for each field of view. The effect appears strong with masking because the variance ctq 
itself becomes smaller for masked maps than for unmasked maps. Unfortunately, it is difficult to 
calculate this effect analytically because the variance depends on a number of elements such as the 
exact geometry and the size distribution of the mask regions, and also on the smoothing scale etc. 
We thus resort to estimating the effect directly using a set of numerical simulations. 

Figure [3] summarizes the effect of mask regions. There, we compare the mean of Vq over 1000 
Gaussian maps with the Gaussian prediction given by Eq. ([6]). For the Gaussian prediction, we 
estimate the quantities (/C), ctq ^^^ ^i from a total of 1000 realizations for masked and unmasked 
maps, respectively. These quantities thus serve as 'global' values. The error bars in each plot 
represent the variance of Vq around the mean. The left and right panels differ in that the MFs 
are plotted as a function of /C — (/C) in the left panel whereas those in the right panel are shown 
as a function of (/C — {IC))/(Tq. Hence the apparent variance of the MFs (error bars) in the right 
panel includes somewhat artificial fluctuations due to the variance of the measured ctq for each 
field of view. In the lower panels, we show the difference between the mean Vq and the Gaussian 
prediction. The difference should be compared with the field variance that are indicated by error 
bars. Notice that the difference from the Gaussian prediction is larger than the field variance when 
the MFs are evaluated with normalization as (/C — (/C))/cro. On the other hand, the Gaussian 
prediction describes well the mean MFs for both masked and unmasked maps as long as the MFs 
are evaluated without normalization of /C by ctq (left panel). Clearly the systematic decrease of o"o 
and 0"! due to mask regions introduce a significant bias to the MFs. 
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Fig. 3. — The effect of normalization of the MFs. We compare the mean Vq over 1000 maps with 
the Gaussian prediction of Eq ([6]). In the right panel, Vq in the left panel is calculated without 
normalization whereas that in the right panel is calculated for each /C field normalized by its variance 
(see text and the horizontal axes label). The gray points in the lower portion show the differences 
between the mean Vq and the Gaussian prediction. The differences are also compared with the 
variance of Vq estimated from our 1000 maps (black error bars), i.e. the standard deviation of Vq 
divided by VlOOO. 



4.2. Ray-tracing lensing simulations 

Having examined the effect of masking on Gaussian fields in the last section, we are now 
able to study the effect in more realistic cases. To this end, we use ray-tracing simulations of weak 
gravitational lensing described in Section [3.21 We pay particular attention to non-Gaussian features 
in the case with masks. The total non-Gaussianity probed by the lensing MFs AV^° is given by 



G{ 



^yobs ^ y. (masked) - F^^ (masked). 



(15) 



where 1^ (masked) is i-th Minkowski Functional on a masked map and VJ*^ (masked) is the Gaussian 
term of V^ (masked). 



We can then decompose AV^° into three components: 
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Fig. 4. — We plot the differences between the lensing MFs on masked ray-tracing simulation maps 
and the Gaussian term. The total non-Gaussianity obtained from the masked maps AV^ (black 
line with closed circle), the non-Gaussianity caused by non- linear gravitational growth AV^^'^'^^' ^ 
(black line with open circle), and the bias of lensing MFs due to mask regions for ray-tracing maps 
^yhias (^gxscy line with closed circle). We also plot the Gaussian term of AT^ '^ (gray line with 
open square). See the definition of each component given by Eq. (fT5|) - p^ . 



where i:\Vf^^^^ ^ represents the non-Gaussianity induced by non-linear gravitational growth, AT^ '^ 
describes the mask bias of MFs for non-Gaussian maps, and AV^ ' corresponds to the Gaussian 
term of AV^'''^^. In order to calculate these quantities, we first need to calculate ^j*^ (masked) and 
VJ (unmasked). For this purpose, we measure the following three quantities from 1000 masked 
ray-tracing maps: 



M = (/C), 



O-n 



(k;2 



^ 



0-1 



IV/CI 



(20) 



The same quantities are measured also for the unmasked lensing maps. We can then estimate 
'V^'^ (masked) and T^j*^ (unmasked) using these quantities and Eq. ([6])- ([8]). For the Gaussian terms. 



we also consider the correction of the finite binning effect pointed out by iLim fc SimonI (J2012l ) . The 
correction is needed because the threshold u to calculate the MFs Vi and V2 is not continuous but 
discrete with some finite width. We calcul ate the correct i on by integrating the analytic formula 
(Eq. ([7]), ([8])) for finite binning width (see iLim &: SimonI (J2012l ) for details). We then estimate 
Vi (masked) and V^ (unmasked) from masked and unmasked maps directly. Figure U] shows the 
various non-Gaussian contributions (Eq. (J16p - ()19p ) calculated directly from 1000 masked ray-tracing 
maps. We find that AV^*^'^^ is comparable to AV^^'^'^^^ ^ in the ray-tracing maps. The mask bias 
^^bias contributes significantly to the observed non-Gaussianity AV^ . Note also that AV^ ^^^' 
is sub-dominant although not negligible for Vi and V2- Clearly the mask bias can be a significant 
contaminant for cosmological parameter estimation using the lensing MFs. In the following, we 
include the bias effect when comparing simulation data and observations. 
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4.3. Application to Subaru Suprime-Cam data 

It is important to test whether we can extract cosmological information from masked noisy 
shear data using the lensing MFs. To this end, we use available Subaru Suprime-Cam data. We 
analyze the observed weak lensing map by using the statistics derived from a large set of ray- 
tracing simulations. We include observational effects directly in our simulations, i.e., mask regions 
and shape noises as described in section [3?2l Figure [5] compares the lensing MFs for the Subaru data 
and those calculated for the ray-tracing simulations. We plot the MFs Vq, Vi, V2 in the top panels. 
In the bottom panels, the thick error bars show the cosmic variance of lensing MFs estimated from 
our 1000 simulated maps, whereas the thin error bars are the sum of the cosmic variance and the 
statistical error. We estimate the statistical error from 1000 randomized realizations, in which the 
ellipticity of each source galaxy is rotated randomly. The statistical error is approximately ~ 1.5 
times the cosmic variance for each bin. In order to quantify the consistency of our results, we 
perform a so-called chi-square analysis. We compute the x^ statistics for the observed lensing MFs, 

^^ = (di - mi)C-\dj - nij) (21) 

where di is the lensing MFs in i-th bin for observation, rrii is the theoretical model, and C is the 
covariance matrix of lensing MFs including the cosmic variance and the statistical error. The cosmic 
variance are estimated from 1000 ray-tracing simulations, and the statistical error are computed 
from 1000 randomized galaxy catalogues. We estimate rrii by averaging the MFs over 1000 ray- 
tracing simulations. We use 10 bins in the range of x = [—3, 3] for each MF. For the binning, we 
have a sufficient number of simulations to estimate the covariance matrix of the lensing MFs. The 
resulting value of x^ per number of freedoms is x^/^dof = 3.35/10,9.69/10,12.8/10 and 29.6/30 
for Vo, Vi, V2 and all the MFs. The analysis includes the cosmic variance and the statistical error 
as well as the effect of mask. We conclude that the observed lensing MFs are consistent with the 
standard ACDM cosmology. 



4.4. Impact of Masking On Cosmological Parameter Estimation 

We study cosmological information content in the lensing MFs with masks. An important 
quantity is the cumulative signal-to- noise ratio S/N for lensing MFs, which is defined by 

{S/Nf = /x*C-V, (22) 

where /^ is a data vector which consists of the lensing MFs Vq, Vi and V2, and C is the covariance 
matrix. In order to calculate (S/N)'^, we construct the data vector from a set of lensing MFs as 

{fii} = {Vo{xi), ..., Vo{xw), Vi(xi), ..., Vi(xio), y2(xi), ..., Falxio)}, (23) 

where Xi{ICi — {lC))/ao is the binned normalized lensing field. 
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Fig. 5. — We compare the observed MFs with those from cosmological ray-tracing simulations. 
In the upper panels, the black line shows the observed lensing MFs and the gray one indicates 
the simulation results. The gray error bars show the cosmic variance obtained from 1000 ray- 
tracing simulations. In the lower panels, the black line shows the difference between the obtained 
MFs and the simulation results. The thin error bars are the sum of the cosmic variance and the 
statistical error while the thick error bars indicate only the cosmic variance. The statistical errors 
are estimated from 1000 randomized galaxy catalogues. 



Figure [6] shows the cumulative signal-to- noise ratio S/N as a function of Xj. Clearly the 
information content is reduced by a factor of two in the case with mask. The degradation is 
explained by the reduced effective area. The solid line shows S/N by scaling C with the effective 
survey area. It closely matches the S/N calculated directly from the masked maps. For a Gaussian 
random field, we expect that the variance of MFs should be inve rsely proportional to the effective 
survey area (e.g. IWinitzki fc KosowskyI Il998l : IWang et al.l l2009l ). The result shown in Figure [6] 
suggests that the effective survey area mainly determines how much cosmological information we 
can gain from weak lensing MFs. 

Let us further quantify the overall impact of bias of the lensing MFs. We perform the following 
simple analysis to investigate the effect of the mask bias on cosmological parameter estimation. For 
each realization r of our simulations, we calculate the x^ value as follows. 



X\r) = i^r) - ^f^-^)C~\^,{r) - ^f -^), 



(24) 



where //i(r) is the estimated lensing MFs from each realization r and fx^ ^°^^ is the theoretical 
template for a given cosmology. In practice, we assume that ^^ ^"'^^ is the average over our 1000 
ray-tracing simulations with or without masks. The lensing MFs ^ii{r) are estimated for each 
masked map, and then we use the covariance matrices of the MFs obtained from a total of 1000 
masked maps. If /ii(r) follows the Gaussian distribution, the distribution of x^('") should follow 
a genuine chi-square distribution. We can then clearly see the impact of bias due to masking on 
cosmological constraints by comparing the resulting distribution of x^('^) for fi- ^°^'^ estimated from 
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Fig. 6. — The cumulative signal-to-noise ratio for the weak lensing MFs. The horizontal axis 
represents the maximum value of binned lensing field used in the calculation of S/N. The open 
circles are the S/N for unmasked 'clean' lensing maps whereas the black points are for masked 
maps. The solid line shows S/N obtained by scaling the covariance matrices of MFs with the 
effective survey area. We adopt the mask regions used for the Subaru Suprime-Cam data (see 
Figure [I]). 



unmasked maps. 

Figure [7] shows the resulting distribution of x^i^) foi' our 1000 masked ray-tracing simulations. 
The black histogram is the probability of X'^(^) for the corresponding model using the average MFs 
over the masked maps whereas the gray one is for the unmasked maps. The thick solid line is 
a genuine chi-square distribution with 30 degrees of freedom, and the dashed line indicates the 
one-sigma region for the x^ values. We find an excellent agreement between the thin histogram and 
the solid line. This means that the binned lensing MFs /Uj(r) can be described well by a Gaussian 
distribution. Interestingly, most of the resulting x^i''") without mask lie outside one sigma regions. 
When we do not take account of bias due to mask regions, 55.3%, 59.4%, 74.9% and 85.4% of the 
realizations lies outside one sigma regions of the chi-square values for Vq, Vi, V2 and all MFs. We 
conclude that the bias of lensing MFs due to mask regions can crucially compromise cosmological 
parameter estimation. 
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Fig. 7. — We compare the distribution of x^(^) foi^ l^i "^™'^ evaluated for 1000 masked maps (black 
histogram) with that for //^ ^°^'^ evaluated for 1000 unmasked maps. The thick solid lines is a 
genuine chi-square distribution with 30 degrees of freedom, and dashed line represents the one 
sigma region. 



5. SUMMARY AND CONCLUSION 

We have used a large number of numerical simulations to examine how mask regions affect 
the lensing MFs by adopting the actual sky-mask used for a Subaru observation. We have then 
compared the observed lensing MFs with the results of cosmological simulations to address whether 
the observed MFs are consistent with the standard cosmological model. 

The weak lensing MFs are affected by the lack of cosmic shear data due mostly to foreground 
contamination. The expected values of the MFs are biased and the variance is increased even for 
a Gaussian convergence field. The bias corresponds to about 5 percent differences of the power 
spectral amplitudes. The bias is induced for the following reasons: (i) masked regions effectively 
reduce the number of sampling Fourier modes of cosmic shear and (ii) masked regions introduce 
scatter of variance of the reconstructed weak lensing mass field for each field of view. The former can 
be corrected analytically as shown in Appendix |Al while numerical simulations are needed to study 
the latter effect. We have used 1000 ray-tracing simulations with mask regions and realistic shape 
noises, to show that the non-Gaussianities detected by the MFs do not solely come from gravity 
induced non-Gaussianities. Masked regions significantly contaminate the pure gravitational signals. 

We have calculated the lensing MFs to the observed weak lensing shear map obtained from 
a Subaru Suprime-Cam imaging survey. Our analysis include, in addition to the cosmic variance, 
the statistical error estimated from 1000 randomized galaxy catalogues. The resulting x^/^dof = 
29.6/30 for all the MFs suggests that the observed MFs are consistent with the adopted standard 
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ACDM cosmology. 

We also perform simple analysis to examine the impact of mask regions on cosmological pa- 
rameter estimation. From the cumulative signal-to-noise ratio for the lensing MFs, we have found 
that the cosmological information content in the MFs can be largely determined by the effective 
survey area. By studying the resulting distribution of x^ value for simulated maps with masks, we 
can characterize how the "mask bias" of the MFs affects cosmological constraints. We have shown 
that most of the resulting x^ values are found outside the expected one sigma region, when the 
mask is not considered. Clearly the mask bias compromise significantly the cosmological parameter 
estimation. 

Finally, we address the ability of the lensing MFs to constrain cosmological models. By as- 
suming a simple scaling of the covariance matrix of MFs by survey area, we can reduce the error 
of MFs in each threshold bin by a factor of ~ 20 (100) for upcoming weak lensing surveys with a 
1000 (20000) deg survey area. In the case of a 1000 deg survey, the error in each bin is translated 
to a ~ 5 percent difference in ug. Similarly, for a LSST-like survey with a 20000 deg^ area, the 
error in each bin is as small as ~ 0.5 percent in erg. The lensing MFs are a promising method for 
cosmology even with mask regions. It is important to model the effect of mask accurately in order 
to make the best use of lensing MFs for cosmological constraints. A simplest way would be to use 
directly the observed mask on ray-tracing simulations as we have done in the present paper. We 
will need such ray-tracing simulations covering a wide area of more than a thousand square degrees 
for lensing MFs in upcoming wider surveys. 

Further extensive studies are needed in order to devise a way to extract pure cosmological infor- 
mation from the lensing MFs. It is important to study ho w much systemati c non-Gaussianities are 



intro duced by, for examp le, source galaxy clustering fe .g. lBernardeaulll998l). source-lens clustering 



(e.g. iHamana et al.ll2002l ). the intrinsic alignment (e.g. iHirata Sz Seljakll2004l ). and inhomogeneous 
ellipiticity noises due to inhomogeneous surface number density of sources (Shirasaki et al., in 
preparation). The upcoming wide-field surveys will provide highly-resolved lensing maps but with 
complicated masked regions. Our study in the present paper may be useful to properly analyze the 
data and to accurately extract cosmological information from them. 
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us with their ray-tracing simulations data. M.S. is supported by Research Fellowships of the Japan 
Society for the Promotion of Science (JSPS) for Young Scientists. This work is supported by World 
Premier International Research Center Initiative (WPI Initiative), MEXT, Japan and in part by 
Grant-in- Aid for Scientific Research from the JSPS Promotion of Science (23540324). Numerical 
computations presented in this paper were in part carried out on the general-purpose PC farm at 
Center for Computational Astrophysics, CfCA, of National Astronomical Observatory of Japan. 
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A. EFFECT OF MASKS ON VARIANCE OF SMOOTHED CONVERGENCE 

FIELD 

Here, we summarize the effect of masked regions on the variance of a smoothed convergence 
field /C. When there are masked regions in a survey area, one needs to follow a special procedure 
in order to construct a smoothed convergence field. Let us define the masked region A^s(0) in a 
survey area as 

, , ,^, I 1 where 6 lies in data region , , , 

I U otherwise. 

When the area with mask A^s(0) is smoothed, there are ill-defined pixels due to the convolution 
between Ms and a filter function for smoothing U{9). We need to discard the ill-defined pixels to 
do statistical analyses. We therefore paste a new mask Mi{6) so that we can mask the ill-defined 
pixels as well. We then get 

/C°^«(6») = 7Wi(6>)/Ci(6>), (A2) 

where 

/Ci(0) = j d^(l)U{e-cj))Ms{cl))K{cj)). (A3) 

The variance of the smoothed field is given by 

i I d^Mm I ^-0^{)C,{£))Cl{£')) exp {^{£ - £') • 6) 

Mii£-£'){ICii£)ICU£')), (A4) 



0"n 



Sj (27r)2(27r)2 
where we use the relation A4i{6)^ = M.i{6). The fourier mode of /Ci is given by 

/Ci(£) = f d'^ed'^^U{e-cl))Msicf))K{cl))exp{i£-e) 

= u{£)J-0^Ms{e'H£-e% 

(/Ci(£)/Ct(0> = U{£)U*{£') J ^^Ms{£i)M:{£[){k{£-£,)k*{£' -£[)) 

= U{£)U*{£') I ^^Ms{£s)M*M){2nf6^'Hi - ii - i' + i[)P.{\i - ii\) 
= U{£)U*{£') j ^Ms{£i)M:{£i+£'-£)P^{\£-£i\) (A5) 
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In the following, we assume that Aii{6) is large enough to cover the ill-defined pixels due 
to smoothing (with a filter function) of the original mask region Ais{6). This means that, with 
Aii{6), there remains only clean regions where the smoothed convergence is not affected by the 
original mask regions A^s(^)- In this case, 

ic"^'{e)^Mi{e) I A^(t)U{e-(i))K{ci)). (A6) 

The fourier mode of /C"^*^ can then be given by 

d'^e Mi{e) I d?cp u{e - 0)^(0) exp (n ■ e) 

d^e'd?ed^(t) Miie')6^^^ {e - G')U{G - (t>)K{cj)) exp {it ■ G) 

dH' 



(2vr)2 

d2f 



d^^'d^^d^,/. Mi{G')U{G - cI))k{cI>) exp {ii ■ G) exp {-il' ■ {G - 6')) 



(2vr)2 

d2f 



Mi{t) I d^ed^cj) U{e - 0)k(0) exp {i{£ - £') ■ G) 
^Mii£')U{£-£')K{£-£'). (A7) 



(27r 
The variance of the smoothed convergence field is calculated as 

\2\ 



4 = \ j d^eiK^'^^e) 

1 /■ d ^ (/C°*^«(£)(/C°^'^)*(£)), (A8) 



Sj (27r)2 
where the ensemble average of the fourier mode is 

(/C°^«(£)(/C°^^)*(£0) ^ j ^^M,{£,)M\{£[)U{£-£,)U*{£' -£',] 

x{k{£-£i)k*{£' -£[)) 
d^h dH[ 



(2vr)2 ^2Tr)^^^^^i)^*^i^i)U{£ - ii)U*{£' - £[) 
x{2TTfP^{\£-£i\)6^^'^{£-£i-£' + £\) 
^Mii£i)Ml{£i +£'- £)\U{£ - £i)\''P^{\£ - £i\). (A9) 

We have checked the validity of Eq. ()A6P by using 1000 Gaussian simulations. They are the 
same set of simulations as in 13.21 For each Gaussian simulation, we paste the observed mask region 
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Fig. 8. — We test the validity of Eq. (|A6p . The gray points with error bars show P {£) obtained 
from 1000 Gaussian maps with mask Ms{6). The gray dashed hne is the theoretical prediction of 
Eq. (IAlip .The black points with error bars show P (i) obtained from 1000 maps without mask. 
The black dashed line is the input power spectrum smoothed by the Gaussian filter U. 



Ais{0) from the Subaru Suprime-Cam observation. The map is then smoothed with a Gaussian 
filter of Eq.Q. The adopted smoothing scale is 1 arcmin. In order to avoid the ill-defined pixels, 
we paste a new mask Aii{0), which is constucted conservatively to cover the regions within 2 times 
the smoothing scale from the boundary of the original mask Ais{0). We then calculate 



pobs^^) 



obs^ 



{jc°'''{i){ic°^y{£))/s. 



If K,"^^ can be well-approximated by Eq. ()A6p . this quantity should be given by 



pobs^^^ ~ 



1 

's 



(2vr) 



UMi{ei)\^\Uii-£i)\^P^{\£-ei\). 



(AlO) 



(All) 



Figure [8] compares Eq. ()A6P and Eq. (jAlip . Clearly Eq. ()A6P is an excellent approximation for the 
observed survey geometry. The ill-defined pixels are efficiently masked by A4i{0). We also find the 
variance ctq decreases by a factor of 0(5%). This causes the bias of Minkowski Functionals even if 
the lensing field is Gaussian. 



